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ABSTRACT 

This report discusses the design and implementation of 
a fixed-based spin simulator and the results derived from 
conducting preliminary spin tests on the simulator. 

The central piece of equipment in the simulator was a 
hybrid computer in which the analog computer solved the 
equations of motion while the digital computer performed the 
tasks of program control and aerodynamic data storage. The 
visual display consisted of a computer-drawn picture on a 
graphics terminal, while pilot control was obtained by use 
of a simulated cockpit situated in front of the graphics 
terminal. 

Results showed that the simulator displayed excellent 
dynamic response characteristics and provided sufficient 
visual cues to perform meaningful spin tests. 

This project was a continuation of previous work and has 
shown that the design and construction of this simulator has 
been an excellent research tool and source for further study 
in the field of control systems and aircraft dynamics. 
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I. 



INTRODUCTION 



The purpose of this project was a continuation study of 
the v;ork done by J. H. Kahrs [Ref. 1] . Kahrs made the 
initial design and construction of a fixed-based variable- 
stability simulator with the task of landing on a carrier or 
runway. The intent here was to make the necessary modifi- 
cations and adjustments to this simulator in order to 
develop a simulator capable of evaluating the spin charac- 
teristics of various aircraft. 

The same basic equipment, i.e. the hybrid computer, 
graphics terminal, and the cockpit simulator, was used as 
constructed by Kahrs. The hybrid's analog computer was used 
for the solution of the equations of motion of the aircraft, 
while the hybrid's digital computer was used for overall 
program control, graphics generation, data reduction and 
storage. The visual display was created by a graphics 
processor and terminal. The cockpit simulator, situated in 
front of the graphics display, consisted of a chair with 
various control levers and buttons v;hose outputs were tied 
directly to the analog computer, thus giving the operator a 
direct means of controlling the simulator. This combination 
of computers and control linkages comprised the basic hard- 
ware of the simulator and formed a complete control loop as 
diagrammed in Figure 1. 
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Figure 1. Simulator Control Loop. 
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The simulator, in the above format, provided the neces- 
sary and required operating and handling characteristics. 
These characteristics, as originally designed by Kahrs , 
included (1) simplicity of setup and operation, (2) ability 
to change aircraft being simulated, (3) capability of simu- 
lating non-linear aerodynamic data, (4) automatic scaling of 
the analog computer solution, (5) fixed-base inside-out 
display, and (6) ability to expand or change simulator task. 

The two most important characteristics desired in this 
project were the ability to expand and change simulator task 
and the ability to change aircraft being simulated. Major 
revisions in the internal programming were needed to account 
for the large differences in the aircraft's flight character- 
istics between the landing phase and the spin phase. These 
differences included the addition of the non-linearized 
equations of motion that are needed in the solution of the 
spin problem. The ultimate goal was to verify that this 
simulator was capable of providing a realistic simulation of 
a spin from a pilot's point of view and at the same time 
obtain meaningful data for further study. 
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II. THE SIMULATOR 



The simulator consisted of three basic pieces of equip- 
ment; the hybrid computer, the graphics terminal, and the 
cockpit simulator. The hybrid computer was used due to the 
outstanding qualities obtained by use of the digital computer 
in conjunction with the analog computer. Due to the rapid 
computational time of the digital computer the bulk of the 
work load was assigned to it. This included the resolution 
of the forces and moments, build-up of graphics data, and 
the non-linearities of the aerodynamic coefficients. The 
analog computer, on the other hand, was used for integration 
of the equations of motion and interfacing. This also 
provided the most rapid means of obtaining a desired output. 
The important feature involved throughout was obtaining the 
desired output data by the fastest means possible. This 
assured that the simulator had an excellent dynamic response 
and led to a display with as little flicker as possible and 
overall graphic quality. 

The graphics terminal consisted of a graphics processor 
and visual display unit (CRT) . Construction of the display 
data was accomplished by the digital computer. These data 
were sent to the graphics processor where a computer-drawn 
picture was formed on the CRT. As new data were compiled by 
the digital computer, the graphics processor would update 
and refresh the visual display, thus providing the pilot 
with a continuous picture. 
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The cockpit simulator was a chair, with attached 
controls, placed in front of the CRT. A Gemini control 
stick was attached to the right arm and provided all three 
control movements: yaw, pitch, and roll. Two programs 

control buttons were provided. One control button was 
located on the top of the stick, while the second control 
button was located on a throttle plate attached to the left 
arm along with a small throttle control. All outputs of the 
cockpit were connected directly to the analog computer for 
processing, thus providing the pilot with the necessary 
controls to run the problem from his station. A complete 
operating manual is contained in Appendix E. 

The simulator was designed .to accept the following input 
data in the form of punched cards: aerodynamic coefficients; 

initial conditions (initial position, altitude, velocity, 
etc.); aircraft constants (weight, moments of inertia, wing 
span, etc.); scale factors based on maximum expected range 
of the problem's variables; and an earth reference grid. 
Appendix F outlines the preparation of this data deck. 

The simulator then automatically processes these data, 
sets the potentiometers of the analog computer, and sets the 
initial conditions on the integrators. Control is then 
transferred to the cockpit where the pilot has three options 
available: 

1. To fly the simulator with the further option of 
aborting the run at any time and returning to fly 
again with the same aircraft parameters. 
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2. To input a new data deck to change the simulated 
aircraft's parameters or to change aircraft being 
simulated. 

3. To stop the complete problem. 

At the completion of each run an analysis of the run is 
briefly displayed on the CRT after which instructions for 
further program control are displayed. The analysis of the 
run consists of the following items: 



1. 


Final 


angle of attack 


2. 


Final 


pitch angle 


3. 


Final 


roll angle 


4. 


Final 


yaw rate 


5. 


Final 


velocity 


6. 


Final 


altitude 
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III. SIX DEGREE OF FREEDOM AIRFRAME EQUATIONS 



Three different coordinate systems — body, wind, and 
stability axes--were available for problem solution. It was 
concluded, on the basis of v;ork done by R. M. Howe [Ref. 2] 
and J. H. Kahrs [Ref. 1], that the best choice of the coor- 
dinate systems was a combination of the wind axes system and 
the body axes system. The translational equations of motion 
would be based on the wind axes system while the rotational 
equations were based on the body axes system. The three 
coordinate systems are shown in Figure 2. The stability and 
body axes system differ by the aircraft's angle of attack, 
while the wind and stability axes system differ by the side- 
slip angle. 

With a knowledge of the three angular rates of motion of 
the aircraft (P,Q,R) about the body axes, normalized values 
for P, Q, and R are obtained. 



P = 



Pb 

2V 



( 1 ) 




( 2 ) 



R = 



Rb 

2V 



(3) 



The angular rates about the stability axes are determined as 
functions of P, Q, and R (Eqn. 4-6). 

P = P cosa + R sina (4) 

s 

Qg = Q (5) 
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Figure 2. Aircraft Coordinate Axes. 
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P sina + R cosa 



( 6 ) 



R = 
s 



The aerodynamic forces expressed in the body axes can 
nov/ be calculated. 
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The body axes forces (Eqn. 7-9) can then be resolved into 
forces along the stability axes.. Thrust (assumed to be in 
the X-direction only) and the force of gravity are also 
inserted . 



X 



m 



T • o . aero 

— - g sin0 + 

m ^ m 



cosa 



g cos0cos({) + 



aero 

m 



sina 



( 10 ) 



y Y 

s Q • A ■ aero 

= g COS0 sin<}) + 



m 



m 



( 11 ) 



m 



T • Q . ^aero 

— - g sin0 + 

m ^ m 



sina 



+ 



g cos0cos4> + 



aero 

m 



cosa 



( 12 ) 



17 



/ 






The final transformation of forces from the stability 



axes to wind 


axes can then be made (Eqn. 13-15) . 
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The translational equations of motion can now be found 
by solving for the derivatives. 



V = - 


— (16) 
in 


-a = 


P sinS 

^ 0 (17) 

mV cosB cosB ^s ^ ^ 


-6 = 


-F 

X 

— 7^ + R (18) 

mV s 



The Euler angular rates are also solved now (Eqn. 19-21) 



-ijj = 


(R cos({) + Q sind) / to \ 

cose 
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-e = 


-Q cos(}) + R sind (20) 
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-P + ijj sin0 (21) 


Next the 


moment equations about the body axes are solved 


for as shov;n 


below . 
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( 23 ) 



JH- = (C + QC + 6itc ) 

I I m m mjc • ^ 

yy yy q <5 it 



(C + PC + RC + 6aC + 6rC ) (24) 

I I n n n ' 

zz zz p r oa or 

The rotational equations of motion can now be written 
(Eqn. 25-27) . 
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(26) 



(I - I ) I (P-QR) 

R = , „yy -. pQ + + N 



zz 



I I 

zz zz 



(27) 



In the development of the velocities in the inertial 
frame, a small angle approximation was not used due to the 
fact that large angles are frequently encountered during the 
spin problem. The resulting equations are as follows: 



S = V 

X 



cosacos6cos0cosip + sing (-coscf)sin>|j + sin(j>sin0cosi/j) 



+ sinacosB (sin<J)siniiJ + cos4>sin0cos>|j) 



(28) 



S = V 

y 



cosacosBcos0sim|^ + sinB (cos(})cosi{)+sinc{)sin0sinip) 



+ sinasinB (sin<|)Cos4^ + cos4>sin0sini|)) 



(29) 



S = V 
z 



-cosacosBsin0 + sinB (sin4>cos0) 
+ sinacosB (cos0cos4>) 



(30) 



Integration of the various derivatives derived in the 
above equations yields the desired state variables. 
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IV. GRAPHICS 



A. GRAPHICS PP^SENTATION 

The need for a good graphics presentation was evident 
due to the complete lack of any physical and audible cues to 
the pilot. Therefore a graphics display had to be generated 
which would give a representation of the earth's surface and 
provide the pilot with enough visual cues to enable him to 
orient himself throughout the problem. The essential visual 
cues needed for the pilot's orientation were indications of 
roll, pitch and yaw. A horizon was included for the orien- 
tation in pitch and roll. Due .to the continuous yawing 
motion of an aircraft during a spin, it was vital to have a 
graphics picture that would accurately display this motion. 
After investigation of several possibilities, it was con- 
cluded that a grid system, attached to the earth's surface, 
composed of a one mile per side squares v;ould provide the 
necessary requirements. The entire grid system was limited 
to six squares per side to avoid excessive computational 
time. By extending the grid lines at all four cardinal 
headings, the visual representation of the ya\<;ing motion of 
the aircraft was further enhanced. This improvement was 
accomplished by the fact that during a continuous yaw the 
pilot only sees a series of lines every 90 degrees. He is 
thus able to positively determine the yaw direction as well 
as a relative measure of the yaw rate. These extensions of 
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the grid lines also assist in the orientation of pitch and 
roll. This grid system is shown in Figure 3. 

To provide the pilot with another source of visual cues, 
a single instrument was constructed. This instrument pro- 
vided the pilot with two sources of information. The first 
was a "needle" providing the pilot with the aircraft's angle 
of bank and the second was a "ball" providing the pilot with 
a measure of the aircraft's sideslip angle. This instrument 
was not the normal "needle-ball" found on all Naval aircraft 
but was a combination of an attitude gyro and turn and bank 
indicator. It was felt that the angle of bank information 
taken from the attitude gyro combined with the sideslip 
angle information taken from the turn and bank indicator 
provided the pilot with the best information in the limited 
space available. 

In the final display, the picture that the pilot saw was 
restricted in direction to the X-axis of the aircraft and to 
a square with a field of view angular limit of plus or minus 
18.5 degrees. Within this square or window lies the horizon 
and portions of the grid reference system. As the aircraft 
maneuvers in space the horizon and grid reference system 
move dynamically within this window creating the sensation 
of flight. The complete initial display is pictured in 
Figure 4 . 

B. GRAPHICS PROCESSING 

The processing of the grid reference system from a three- 
dimensional object to a two-dimensional plane capable of 
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Figure 4. Initial Graphics Display. 
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being displayed on a screen proved to be the most complicated 
portion of the display. The same system developed by 
L. G. Roberts [Ref. 3] and R. B. Desens [Ref. 4] and used by 
J. H. Kahrs [Ref. 1] was used here. It basically consists 
of the construction of a single transformation matrix for 
all display points. Two coordinate systems were used in the 
development of this matrix, that of the earth reference 
system (object) and that of the aircraft (viewing plane) , 
and are shown in Figure 5 . 

The transformation matrix (H-Matrix) is the product of 
five matrices: rotation, translation, perspective, scale, 

and another translation. The H-Matrix is further reduced 




Figure 5. Graphic Coordinate System. 



by Desens to the product of three matrices as shown in 
Figure 6. 
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Figure 6. H-Matrix. 



The first matrix consists of the direction cosines of 
the Euler angle rotation of the object axes and the trans- 
lation distances between the two axes system. The direction 
cosines and the order of rotation are shown in Figure 7. 

Since the object axes in this problem was a fixed earth 
reference system, the direction cosines are all zero except 
for the diagonal terms which have a value of one. The 
matrix then is simply reduced to a viewing plane translation. 
The translation distances; X^p, Y^ and are derived as 
follows : 

^T ^object ^view plane 

Y = Y - Y 

T object view plane 

Z r= Z — Z 

T object view plane 
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A 2 = cos3sina 
= - sin3 
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B 2 = cosacosy + sinysinBsina 
B^ = sinycosB 

= sinysina + cosysingcosa 
C 2 = - sinycosa + cosysingsina 
= cosgcosy 

Order of rotation a,3/y 
Figure 7. Euler Angle Rotation. 

The second matrix is composed of the direction cosines 
of the Euler angles of the aircraft or viewing plane axes 
and are listed below. 



AAj^ = 


cosi|;Cos6 
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CC 2 = - cos0sin({) 

CC^ = cos0cos<f) 

Order of rotation 

ip (yaw) , 0 (pitch) , (t>(roll) 

The third matrix is composed of the offset option, 
and , the scale factor (S) and the focal length (F) . The 
offset option was not used but its usage is explained in 
Ref. 4. The scale factor was set at 1/2 and the focal 
length, the distance between viewer and the viewing plane, 
was set at 1.5 ft. 

To process the grid reference system or object, a point 
(X,Y,Z) is taken and converted to homogeneous coordinates by 
the addition of a scale factor (W) . This yields new coor- 
dinates (X',Y',Z',W) where X' = WX, Y' = V7Y , Z' = WZ . Post 
multiplying these coordinates by the H-Matrix yields 
(X" , Y" , Z" ,W ) . The display coordinates can now be found by 
dividing by W . 

Y = Y"/W 

Z = Z"/W 

The X-coordinate was not used in this display since it 
is only an indication of depth. The final display coordin- 
ates are obtained by passing each point of the grid reference 
system through the H-Matrix. 

The horizon was constructed by drawing a line parallel 
to the Y-object axis and offset a distance in the negative 
X-direction. The Z-coordinate was set to zero. In doing so 
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the assumption was made that the earth was flat and there- 
fore did not account for the fact that when flying at an 
altitude a depression angle is created due to the curvature 
of the earth. This assumption was valid since the horizon 
was only used for a relative measure of orientation for the 
pilot and not as a measure of the altitude of the aircraft. 
The coordinates of the horizon are passed through a modified 
viewing plane orientation matrix to account for pitch and 
roll of the aircraft. The yaw angle of the aircraft is set 
to zero since it does not have any apparent visual influence 
on the movement of the horizon. This results in the follow- 
ing matrix. 

AAj^ = COS0 

AA 2 = - sin0sin(}) 

AA^ = sin0cos<}) 

BB^ = 0 
BB 2 = coscf) 

BB^ = sincf) 

CC^ = - sin0 
CC 2 - cos0sin(j) 

CC^ = cos0cos(J) 
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V. RESULTS OF SPIN TESTS 



Several spin tests were conducted in order to compare 
these results with the results of other simulators using the 
same aircraft parameters and aerodynamic coefficients. The 
aircraft parameters and aerodynamic coefficients used in all 
tests were obtained from W. P. Gilbert [Ref. 5] and are 
representative of a variable-sweep fighter aircraft. A wing- 
sweep angle of 16° was the only configuration considered 
during the tests. 

In all cases the aircraft was initially positioned at an 
altitude of 30/000 feet and an .airspeed of 622 ft/sec. The 
spin was entered by reducing the throttle to the idle posi- 
tion followed by movement of the horizontal stabilator to 
the full trailing edge up, rudder to full trailing edge left, 
and the ailerons to full right wing down. The controls were 
maintained in these positions until the spin had fully 
developed when they were then released to the neutral posi- 
tion. The resulting spin was limited in altitude to 10,000 
feet or until the spin development was such that the limits 
of the computer were exceeded. At no time was an attempt 
made to effect a recovery. 

The resulting aircraft motion was a stall exhibiting 
unstable flight characteristics in lateral control. A right 
or left roll of 360° at high angles of attack then ensued 
after which the aircraft entered a fast flat-spin. During 
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the spin, the average angle of attack was 80° with a yaw 
rate of about -133 deg/sec. The pitch angle and roll angle 
remained relatively constant throughout the spin at approx- 
imately -15° and 0° respectively. At the end of 60 seconds, 
the aircraft had completed 9.5 turns and lost about 7000 
feet of altitude. The final aircraft velocity was 262 
ft/sec. A complete time history of the results of one of 
the spin tests is shown in Figure 8. A brief comparison 
with the results obtained by Gilbert is shown below. A 
complete time history of Gilbert's results is shown in 
Figure 9. 

Simulator Gilbert 
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Figure 8. Results of Spin Test. 
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Figure 8. (Continued). 
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Figure 8. (Continued). 
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Representative Spin Results 



VI. CONCLUSIONS 



As a result of conducting a series of spin tests on this 
simulator the following conclusions have been reached: 

1. The simulator provided excellent visual cues to 
enable experienced pilots to orient themselves 
throughout the entire problem. 

2. The simulator exhibited good dynamic response with a 
computational loop time of 80 milliseconds and solu- 

Z' 

tion update rate of approximately 16 samples per 
second. 

3. The simulator v/as capable of providing meaningful 
spin data from which comparisons and further study 
could be made. 

4. One minor drawback to the simulator was the lack of 
any physical cues, i.e. "g" forces. These forces can 
be significant in an actual spin and severely res- 
trict the pilot's movement. Lack of these forces 
has led to pilots over-controlling the simulator 

and thereby risking the possibility of introducing 
erroneous information or exceeding the limits of the 
computer . 

5. The versatility of the simulator, in the ease with 
which the aircraft or aircraft parameters may be 
changed, has made it a valuable research tool and 
source for further study in the field of flight 
dynamics . 
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APPENDIX A 



THE DIGITAL PROGRAM 

This appendix contains a listing of the computer program 
written in the FORTRAN language. It consists of a main 
program and three subprograms: the data reduction, the 

computation loop, and the A/D-D/A. The A/D-D/A subroutine 
is written in METASYMBOL, the assembly language for the 
XDS-9300, and was' -prepared by the computer laboratory staff. 
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APPENDIX B 



THE ANALOG PROGRAM 
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APPENDIX C 



PROGRAM VARIABLES 

The following is an alphabetical listing of the program 
variables used in the computer program. 



A 


Array of angles, in order, alpha, beta, psi, theta 
and phi 


AlO 


Initial conditions on A(l), the angle of attack 


AA 


Array of digital to analog variables 


AD 


Array of analog to digital variables 


ADOTN 


Array of the negatives of the time derivatives of 
the angles in A 


ADO TO 


Old derivatives of angles in A used for prediction 


ALI 


Aerodynamic rolling moment about the body axis 


ALPHA 


Real argument based on angle of attack used for 
coefficient lookup 


AMI 


Aerodynamic pitching moment about the body axis 


AN I 


Aerodynamic yawing moment about the body axis 


ARG 


Argument used during sine and cosine lookup (type 
real) 


B 


Wing span 


B2 


One half the wing span B 


BETA 


Real argument based on sideslip angle used in 
coefficient lookup 


BX 


A measure of sideslip angle used in the generation 
of the "ball" 


C 


Array of aerodynamic coefficients for a given angle 
of attack and sideslip angle 


CA 


Array of cosines of the angle array A 
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CB 


Mean aerodynamic chord 


CB2 


One half the mean aerodynamic chord CB 


CCA(n) 


Scalar equivalents of the array CA 


CHOYC 


Choice of program options offered by main program 


COEFA 


Array of coefficients dependent on angle of attack 
only 


COEFAB 


Array of coefficients dependent on angle of attack 
and sideslip angle 


COF (n) 


Summation variables used in the interpolation 
process of coefficient lookup 


CONI 


Program constant equal to the number of degrees in 
a radian . 


DA 


Aileron deflection 


DEGA 


Degrees of angle of attack for which COEFA and 
COEFAB are tabulated 


DEGB 


Degrees of sideslip for which COEFAB are tabulated 


DIT 


Flying tail deflection 


DITO 


Initial flying tail deflection setting 


DIV 


Scale factor divisor used to divide out the scale 
factor of the homogeneous transformation coordinates 


DR 


Rudder deflection 


F 


Focal length of viewing image 


FXAM 


Aerodynamic force in the X-direction divided by 
the mass of the aircraft 


FXSM 


Force in X-direction of the stability axis divided 
by the mass of the aircraft 


FYAM 


Aerodynamic force in the Y-direction divided by 
the mass of the aircraft 


FYSM 


Force in the Y-direction of the stability axis 
divided by the mass of the aircraft 


FYWM 


Force in the Y-direction of the wind axis divided 
by the mass of the aircraft 



,80 



f- 




'■4 

• a 



• j 



. 'I 






FZAM 


Aerodynamic force in the Z-direction divided by 
the mass of the aircraft 


FZWM 


Force in the Z-direction of the wind axis divided 
by the mass of the aircraft 


G 


Acceleration due to gravity 


GEE 


Scaled acceleration due to gravity 


HINSCT 


Horizontal intersection of a line, used in the 
software window 


HM 


Viewing plane orientation matrix for the horizon 


I 


Integer variable usually a counter 


lA 


Do-loop counter based on angle of attack 


lAGN 


A counter to check the number of times a point has 
gone through the software window 


I ALP HA 


Integer conversion of ALPHA 


lARG 


Integer conversion of ARG 


IB 


Do-loop counter based on the sideslip angle 


I BETA 


Integer conversion of BETA 


ICHOYC 


Integer conversion of CHOYC 


IDE 


Graphic array containing the dynamic portion of 
the display 


I DEV 


Graphic device number (1 or 2) 


IDIR 


Graphic block directory for the graphics digital 
processor 


lER 


Error flag returned by graphic subroutines 


IPLS 


Integer counter equal to I plus one 


ISF7 


Integer scale factor based on thrust divided by 
mass 


ISQ 


Array/graphics block of fixed data, basically the 
square or window 


ITDIR 


Text directory for the graphics digital processor 


ITEXnn 


Line of text that is sent from the digital computer 
to the graphic digital processor 



81 



IX 

IXl 

1X2 

lY 

lYl 

IY2 

J 

K 

KK 

KPOT 

LTR 

NAD 

NDA 

NULL 

NUM 

P 

PB 

PDOTN 

POT 

PS 

Q 

QB 

QDOTN 

QS 

QUE 



Operating values of start/end points (X-coordinate) 
used in software window 

Equivalent of XSTART, used in software window 

Equivalent of XEND , used in software window 

Operating values of start/end points (Y-coordinate) 
used in software window 

Equivalent of YSTART, used in software window 
Equivalent of YEND, used in software window 
Do-loop counter 
Do-loop counter 

Array of integers used in outputing coefficients 

Array of integers used in outputing analog pot 
settings 

Clock variable used in loop timing 

Number of analog to digital conversion variables/ 
truck lines 

Number of digital to analog conversion variables/ 
truck lines 

Octal value used to null out a line of text 
Variable used for program flow and control 
Angular velocity of the aircraft about X-axis 
P normalized 

The negative of the time derivative of P 
Array of pot settings for the analog 
P about the stability axis 

Angular velocity of the aircraft about the Y-axis 
Q normalized 

The negative of the time derivative of Q 
Q about the stability axis 
Dynamic pressure 
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R 

RB 

RDOTN 

RHO 

RIXX 

RIXZ 

RIYY 

RIZZ 

RMASS 

RS 

RT 

S 

SA 

SDOT 

SF 

SFl 

SF2 

SF3 

SF4 

SF5 

SF6 

SF7 

SF8 



Angular velocity about the Z-axis 
R normalized 

The negative of the time derivative of R 
Density of air 

Moment of inertia of the aircraft about the X-axis 

Product of inertia about the X and Z-axis 

Moment of inertia of the aircraft about the Y-axis 

Moment of inertia of the aircraft about the Z-axis 

Mass of the aircraft 

R about the stability axis 

Graphic transformation matrix or array 

Wing area of the aircraft 

Array of sines of the angle array A 

Array of time derivatives (velocities) of the 
distances in all three inertial directions 

Product of the scale factor and focal length 

Scale factor based on maximum velocity in 
X-direction 

Scale factor based on maximum velocity in 
Y-direction 

Scale factor based on maximum velocity in 
Z-direction 

Scale factor based on maximum angular velocity in 
X-direction 

Scale factor based on maximum angular velocity in 
Y-direction 

Scale factor based on the maximum angular velocity 
in the Z-direction 

Scale factor based on thrust divided by mass of the 
aircraft 

Scale factor based on maximum angle of attack 
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SF9 


Scale 


factor 


based 


on 


the 


maximum 


sideslip angle 


SFIO 


Scale 


factor 


based 


on 


the 


maximum 


yaw angle 




SFll 


Scale 


factor. 


based 


on 


the 


maximum 


pitch angle 




SF12 


Scale 


factor 


based 


on 


the 


maximum 


roll angle 




SF13 


Scale factor 
X-direction 


based 


on 


the 


maximum 


distance in 


the 


SF14 


Scale factor 
Y-direction 


based 


on 


the 


maximum 


distance in 


the 


SF15 


Scale factor 
Z-direction 


based 


on 


the 


maximum 


distance in 


the 


SF16 


Scale factor based 
deflection angle 


on 


the 


maximum 


flying tail 




SF17 


Scale factor 
tion angle 


based 


on 


the 


maximum 


aileron deflec- 


SF18 


Scale factor 
tion angle 


based 


on 


the 


maximum 


rudder deflec- 


SF19 


Scale 


factor 


based 


on 


the 


density 


of air 




SF20 


Scale 


factor 


based 


on 


SFl 


squared 


divided by 


SF19 


SF21 


Ratio 


of SF6 


to SF4 












SF22 


Ratio 


of SFl 


to SF4 










SF23 


Ratio 


of SFl 


to SF5 










SF2 4 


Ratio 


of SFl 


to SF6 










SF25 


Ratio 


of SF4 


to SF5 










SF26 


Ratio 


of SF5 


to SF6 










SF27 


SF4 squared 














SF28 


SF5 squared 














SF29 


SF6 squared 














SF30 


Feedback gain for yaw 


damper 






SF31 


Scale 


factor 


based 


on 


SF7 


divided 


by SFl and 


SF5 


SF32 


Scale 


factor 


based 


on 


SF7 


divided 


by SFl and 


SF6 
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SF33 


Ratio of SFl to SF2 


SF34 


Ratio of SFl to SF3 


SF35 


Scale factor based on ratio of SF5 to SF8 


SF36 


Scale factor based on ratio of SF6 to SF9 


SF37 


Scale factor based on ratio of SF6 to SFIO 


SF38 


Scale factor based on ratio of SF6 to SFll 


SF39 


Scale factor based on ratio of SF6 to SF12 


SF40 


Ratio of SF7 to SFl 


SFA 


Array of scale factors, made up of SF8 thru SF12 



SFACTOR Scale factor of the H-Matrix 



SK 


Variable indicating portion of each line to be 
removed from display in the window chop loop 


SLOPE 


Slope of a display line, used in software window 


SN 


Negative of the scale factor 


SSAn 


Scalar equivalents of SA 


SX 


Actual distance of aircraft in X-direction 
(inertial) 


SXO 


Initial condition on SX 


SXN 


Normalized distance of aircraft in X-direction 
(inertial) 


SY 


Actual distance of aircraft in Y-direction 
(inertial) 


SYN 


Normalized distance of aircraft in Y-direction 
(inertial) 


SZ 


Actual distance of aircraft in Z-direction 
(inertial) 


SZO 


Initial condition on SZ 


SZF 


Final cutoff value on SZ (10,000 ft.) 


SZN 


Normalized distance of aircraft in Z-direction 
(inertial) 


T 


Thrust of the aircraft 
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TAU 


Time delay for the prediction and update 


TLINE 


Array of the earth grid reference points 


TM 


Thrust divided by mass of aircraft 


TRIG 


Array of sines of various angles 


TRIGC 


Array of cosines of various angles 


TX 


A measure of roll angle in the X-direction used 
in generation of the "needle" 


TY 


A measure of roll angle in the Y-direction used 
in generation of the "needle" 


V 


Velocity 


VO 


Initial condition on velocity V 


VDOT 


Time derivative of V 


VDOTO 


VDOT of the iteration before 


VINSCT 


Vertical intersection of a line used in window 


VX 


Velocity in the X-direction 


VY 


Velocity in the Y-direction 


VZ 


Velocity in the Z-direction 


W 


Weight of the aircraft 


X 


X value used in software window 


XEND 


Array of X-coordinate ending points of display lines 


X START 


Array of X-coordinate starting points of display 
lines 


XTE 


X-coordinate ending point for the horizon 


XTEMP 


Temporary storage for X 


XTS 


X-coordinate starting point for the horizon 


Y 


Y value used in software window 


YEND 


Array of Y-coordinate ending points of display lines 


YSTART 


Array of Y-coordinate starting points of display 
lines 
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YTE 


Y-coordinate ending point for the horizon 


YTEMP 


Temporary storage for Y-coordinate in softvvare 
window 


YTS 


Y-coordinate starting point for the horizon 


ZEND 


Array of Z-coordinate ending points of display lines 


Z START 


Array of Z-coordinate starting points of display 
lines 


ZTE 


Z-coordinate ending point for the horizon 


ZTS 


Z-coordinate starting point for the horizon 
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APPENDIX D 



DATA DECK PREPARATION 

The following is a guide for use in preparing a data 
deck for use with the simulator. 

Cards 1-3 - These cards contain the inputed scale 
factors as defined in Appendix C. Each scale factor has a 
field width of ten and are ordered as follows: 

Card 1 - SFl, SF2 , SF3 , SF4 , SF5, SF6, SF8, SF9 

Card 2 - SFIO, SFll, SF12 , SF13, SF14, SF15, SF16, SE17 

Card 3 - SF18, SF19 , SF30 , the remaining columns are 
blank 

Card 4 - This card contains the initial conditions of 
the problem. Each number has a field width of ten and they 
are ordered as follows: initial X-coordinate (SXO) in feet, 

initial altitude (SZO) in feet, initial velocity (VO) in 
ft/sec, initial angle of attack (AlO) in degrees, initial 
elevator angle (DITO) in degrees and the density of the air 

3 

(RHO) in slug/ft . The remaining portion of the card is 
blank. 

Cards 5 and 6 - These cards contain the aircraft con- 
stants. Each constant has a field width of ten and they are 
ordered as follows: 

Card 5 - Weight (W) of aircraft in lb. , wing span (B) in 

ft. , wing chord (CB) in ft. , wing area (S) in 

ft. , moment of inertia about X-axis (RIXX) in 
2 

slug-ft , moment of inertia about the Y-axis 
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2 

(RIYY) in slug-ft , moment of inertia about the 
. . 2 

Z-axis (RIZZ) in slug-ft , and the cross product 
of inertia (RIXZ) in slug-ft^. 

Card 6 - Maximum thrust (T) in lbs . , altitude cutoff 

(SZF) in feet. The remaining columns are left 



blank . 

Cards 7 - 412 - These cards contain the aerodynamic 
coefficients dependent on both alpha and beta. The first 
card in each of a sequence of 135 three card series is as 
follows: in columns 1 and 2 the coefficient number, columns 

3 and 4 blank, column 5 the index on beta (IB) , columns 6-10 
blank. The remaining portion of the card is filled with 
coefficients indexed on alpha (lA) in ten column intervals. 
The second and third cards in the series contain the remain- 
ing coefficients in ten column intervals and indexed on 
alpha. Alpha is indexed 19 times for each beta index while 
beta is indexed 9 times for each of the fifteen coefficients. 
The range on alpha is 0 - 90 degrees with 5 degree intervals 
while the range on beta is -40 degrees to +40 degrees with 
10 degree intervals. The coefficients are as follows: 



Number 

1 

2 

3 

4 

5 



Coefficient 
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Number 



Coefficient 



6 

7 

8 
9 

la 

11 

12 

13 

14 

15 
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Cards 413 - 431 - These cards contain the aerodynamic 
coefficients dependent only on alpha. The first card in each 
of a sequence of 6 three card series is as follows: columns 

1 and 2 the coefficient number, columns 3-10 blank, columns 
11-80 the coefficients indexed on alpha and in ten column 
intervals. The second and third cards contain the remaining 
coefficients again in ten column intervals and indexed on 
alpha. The coefficients are as follows: 



Number 


Coef f i 
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Cards 432-440 - These cards also contain aerodynamic 
coefficients based on alpha only. The only difference in 
the setup of these cards in the field width is increased to 
eleven. As a result the first card in the series is as 
follows; columns 1 and 2 the coefficient number, column 3 
blank, columns 4-80 contain the coefficients indexed on 
alpha and in eleven column intervals. The second and third 
cards contain the remaining coefficients indexed on alpha 
and in eleven column intervals. The coefficients are as 
follows: 



Number 

22 

23 

24 



Coefficient 
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Cards 440 -454 - These cards contain the actual coordi- 
nates for each line of the earth grid reference system. 

Each card represents one line with its starting and ending 
coordinates. On each card the first eight columns are the 
starting X-coordinate , columns 9-16 the starting Y-coordi- 
nate, columns 17-24 the starting Z-coordinate , columns 25-32 
the homogeneous coordinate W, columns 33-40 the ending 
X-coordinate, columns 41-48 the ending Y-coordinate , columns 
49-56 the ending Z-coordinate, and columns 57-64 the homo- 
geneous coordinate W. The remainder of the card is left 
blank. 
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APPENDIX E 



OPERATING MANUAL 

T\i7o tapes are available for use in execution of the 
digital program. One tape (labeled SPIN SIMULATOR SI) 
requires approximately 15 minutes to compile prior to execu- 
tion. The second tape (labeled SPIN SIMULATOR CD) , which is 
a core dump of the first tape, requires only 10 seconds prior 
to execution. It .is recommended that the second tape be used 
due to its short compilation time unless a printed output of 
the program is desired. 

All of the required cards, )30th control and data, are 
located in the laboratory files under FIXED-BASE SPIN 
SIMULATOR. Also filed here is the master card version of 
the program and should not be used unless the tapes are 
destroyed. 

Detailed operating instructions for both the hybrid and 
graphics computers are available in the laboratory office. 

The following instructions should be followed in order to 
insure that the simulator is set up and executing properly. 

1. Load the tape selected on the tape drives. Select 
the proper mount (2 for SI, 3 for CD). 

2. Load and ready the card reader with control cards 
followed by the DATA. 

3. Ready the Line Printer. 

4. On the Digital Main Frame, Punch IDLE, RESET, RUN, 
CARDS . 
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5. Load the Analog Program board and the Analog Logic 
board into the computer (Boards number 7) . 

6. On the Analog Keyboard punch KEYBOARD, LOCAL, POTSET 

In the far right control box, only the following can 
be lit: IDCXl or IDCS.l and REAL TIME. All other 

lights must be punched out. When complete punch 
DIGITAL CMPTR. 

7. On the Analog the following options exist: 

DSO - UP - no printed output 

CENTER - printed output 
DSl - UP - distance integrators disabled 
CENTER - distance integrators enabled 
DS2 - UP - yaw damper ON 

CENTER - yaw damper OFF 

DOWN positions on these switches are mementary con- 
tact positions. 

8. At the Graphics terminal teletype to be used, load 
the controling program by typing 

RESET ("GATDl", 104) ! 

Then execute the program by typing GATED! At this 
point nothing should be on the display screen. 

9. By this time the Digital computer should have 
compiled and loaded the main program and the follov/- 
ing will be typed out on the digital consol teletype 

INPUT AGT NUMBER 

The input light will then come on. To respond, type 
the AGT number (1 or 2) of the unit being used, 
followed by a carriage return. 
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10 . 



The Title should appear on the selected AGT. 

11. The Data will be read in from the card reader. 

12. If output is selected, the output will be printed 
on the line printer. 

13. The Analog Pots will be set by the Digital computer. 
The Address and Ratiometer readings should be 
changing. 

14. Set the cockpit in front of the AGT (will only 
reach to AGT-2). The cockpit is stored behind the 
Analog computer. 

15. The Instructions should flash on the screen. 

16. Punch the button on the Throttle plate and the 
display will appear. 

17. Punch the button on the control stick to fly. If 
during a run you wish to abort that run, punch the 
button on the Throttle plate. 

18. At the completion of the run the Spin Results will 
be displayed, followed by a short delay in which to 
read them. 

19. At this time two options are offered. To fly again 
punch the button on the Control stick and the pro- 
gram jumps back to the Instructions (#15) . 

20. To receive the expanded program options, punch the 
Throttle button. Follow the displayed instructions 
and type your selection on the AGT teletype unit. 
This will result in one of three things happening. 
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1. Logical STOP to the program. 

2. Jump to the Instructions (#15). 

3. Jump to the selection of the ACT number (#9). 
21. At the completion of the computer time, return all 

used equipment (cards, tape, boards, cockpit) to 
their storage locations. 
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APPENDIX F 



SAMPLE OUTPUT 

The following is a sample output of the digital program 
which contains the following sections: 

1. Scale Factors and Aircraft Constants 

2. Original Aircraft Aerodynamic Coefficients 

3. Usable Aircraft Aerodynamic Coefficients 

4. Pot Settings for the Analog Computer 

5. Output of the Amplifiers and D/A Trunks 

6. Earth Grid Reference Lines 
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